FEAT: expm1
ufunc implementation
#183
Merged
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Copilot Summary
This pull request adds support for the
expm1
function (which computesexp(x) - 1
with high precision) to the quad-precision and long double numerical types in the project. The implementation includes the core operation, integration with the NumPy interface, updates to documentation, and comprehensive tests to ensure correctness and precision.Functionality additions:
quad_expm1
andld_expm1
functions to computeexpm1
for quad-precision (Sleef_quad
) and long double types inops.hpp
. [1] [2]expm1
function with the NumPy ufunc interface, making it available for use in Python via the quad-precision type.Documentation and testing:
expm1
is now implemented and tested for both quad and long double types.expm1
intest_quaddtype.py
, covering a wide range of values including edge cases, small inputs, large inputs, and special values (inf, nan).